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C " 3 ' We review recent results based on an application of the real-space renormalization group 

l/^ (RG) approach to a network model for the integer quantum Hall (QH) transition. We 

C D , demonstrate that this RG approach reproduces the critical distribution of the power 

transmission coefficients, i.e., two-terminal conductances, Pc(G), with very high accuracy. 
' The RG flow of P(G) at energies away from the transition yields a value of the critical 

exponent v that agrees with most accurate large-size lattice simulations. A description of 
how to obtain other relevant transport coefficients such as i?L and i?H is given. From the 
1 non-trivial fixed point of the RG flow we extract the critical level-spacing distribution 

' (LSD). This distribution is close, but distinctively different from the earlier large-scale 

Q ' simulations. We find that the LSD obeys scaling behavior around the QH transition with 

^ , V = 2.37±0.02. Away from the transition it crosses over towards the Poisson distribution. 

We next investigate the plateau-to-insulator transition at strong magnetic field. For a 
fully quantum coherent situation, we find a quantized Hall insulator with Rn f» h/e^ 
' up to i?L ~ 20/i/e^ when interpreting the results in terms of most probable value of 

I^N , the distribution function P{Rii). Upon further increasing i?L — * oo, the Hall insulator 

^ ■ with diverging Hall resistance ijjj oc ij£ is seen. The crossover between these two regimes 

depends on the precise nature of the averaging procedure for the distributions P(Ri^) and 
P(-Rh)- We also study the effect of long-ranged inhomogeneities on the critical properties 
of the QH transition. Inhomogeneities are modeled by a smooth random potential with 
a correlator which falls off with distance as a power law r~°'. Similar to the classical 
percolation, we observe an enhancement of u with decreasing a. These results exemplify 
the surprising fact that a small RG unit, containing only five nodes, accurately captures 
most of the correlations responsible for the localization-delocalization transition. 

Keywords: integer quantum Hall effect; real-space renormalization group approach; net- 
work models; energy-level statistics; macroscopic inhomogeneities; quantum Hall insula- 
tor. 
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1. Introduction 

Measuring the resistance of a two-dimensional (2D) electron gas at very low tem- 
perature and subject to a strong perpendicular magnetic field B reveals a striking 
macroscopic quantum phenomenon [1]: the Hall resistance i?H shows very precise 
plateaus at {l/N)h/e'^ where N is an integer number. In the same region the vanish- 
ing longitudinal resistance i?L indicates a dissipationless flow of current. Quantiza- 
tion can also be observed at fractional values of N [2] , which leads to the distinction 
between the fractional and the integer quantum Hall (QH) effect. Despite the simi- 
larities regarding the experimental observations the theoretical description of both 
effects differs considerably. The integer QH effect can be explained reasonably well 
within a non-interacting electron picture, while interactions play a fundamental role 
in the fractional QH effect [3,4]. 

The integer QH effect can be rationalized at a phenomenological level as a se- 
ries of localization-delocalization transitions [5]. For non- interacting electrons, the 
density of states consists of distinct Landau bands. These are broadened by the 
intrinsic disorder of the QH samples. Each such band contains localized states in 
the tails and extended states in the center of the band. As the Fermi energy passes 
through the extended states, the Hall resistance i?H rises to the next plateau and 
the longitudinal resistance i?L is non- vanishing. In the range between the bands 
and for the localized states the Hall resistance remains fixed at the plateau value 
as shown in Fig. 1. 
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Fig. 1. Schematic behavior of longitudinal {axx) and Hall {cFxy) conductivities as well as the 
density-of-states (DOS) for the integer QH effect. _Ep denotes the Fermi energy and the energy of 
the Landau levels is E = hjj^^N — 1/2), N = 1,2, . . lo^ = eB/m the cyclotron frequency. The 
width of the region of extended states at the band centers of the disorder-broadened Landau levels 
is finite due to finite temperatures and finite system sizes. 
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Several theoretical approaches have been proposed to explain the sequence of 
integer QH transitions in more detail [3-10]. For example, it has been established 

that the QH transition is a second-order phase transition [11] where a power-law 
divergence of characteristic length scales is observed. The divergence, e.g. of the 
localization length A of the wave functions, A oc €~'^, where e corresponds to the 
distance from the transition, can be quantified by the exponent ly. Its value is a 
universal quantity, independent from microscopic details of electron motion, disorder 
realization and even the chosen material. 

The Chalker-Coddington (CC) model is a convenient representation of the QH 
situation in the fully quantum coherent regime at low temperature and strong mag- 
netic field [12]. The model describes a single Landau band and thus contains only 
one QH transition. Electrons are treated semiclassically as non-interacting particles 
moving in a smoothly varying 2D disorder potential under a strong perpendicular 
magnetic field. The percolation of the electron through the sample is then given by 
motion along equipotentials interrupted by scattering events at saddle points (SPs) 
of the potential. For theoretical investigations, equipotentials and SPs are mapped 
on links and nodes of a 2D network, respectively, and large networks have to be 
studied in order to reduce finite-size effects [12-17]. 

In the present review, we concentrate on an alternative approach to the CC 
model. Namely, we shall introduce a real-space renormalization-group (RG) ap- 
proach [18,19], based on mapping a characteristic part of the CC network — the RG 
unit — consisting of several SPs onto a single super-SP residing in a new CC-super 
network. This approach is motivated by the success of RG techniques in statistical 
physics and particularly the study of phase transitions [20,21]. Common to both 
real-space [22] and field-theoretical [23] RGs is the elimination of irrelevant (either 
short-range or short-period) degrees of freedom. Afterwards the original phase vol- 
ume is restored by a scale transformation. The system has now the same structure 
as the original but the parameters, e.g. coupling constants, are renormalized. During 
the iteration of the above RG steps the parameters approach their fixed point (FP) 
values and an analysis of this convergence allows one to derive the critical properties 
of the transition. Large effective system sizes can be achieved in this way for the 
CC model. This advantage is counterbalanced by the approximate description of 
the network due to the limited size of the RG unit. 

The RG procedure provides a natural route to studying the complete distribution 
functions of the 2D resistivities defined via 



where Pxx is called longitudinal and pxy Hall resistivity. The conductivity tensor is 
defined as the inverse of the resistivity tensor which yields 
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We emphasize that in 2D resistances (conductances) and resistivities (conductivi- 
ties) are related by a simple, dimensionless geometrical factor. 

As a first application of the real-space RG approach to the CC model, wc study 
the distribution functions P of two-point conductance G and the resistances -Rl 
and Ru- From these, we compute the value of the critical exponent v. Furthermore, 
closing the incoming and outgoing links of the RG unit and thus quantizing the 
allowed energies, wc can make contact with previous stiidics of energy level statistics. 
A suitable scaling ansatz again allows the independent estimation of v. These results 
are in good agreement with previous studies and validate the present RG approach. 

A further application is concerned with the possible existence of the quantized 
Hall insulator. While it is established that plateau-plateau and insulator-plateau 
transitions exhibit the same critical behavior [24-34] the value of the Hall resistance 
i?H in this insulating phase is still rather controversial [35]. Various experiments 
have found that i?H remains very close to its quantized value /i/e^ even deep in 
the insulating regime when already Rl ^ h/e'^ [25-28,36]. This scenario has been 
dubbed the quantized Hall insulator. On the other hand, theoretical predictions 
based on quantum coherent models show that a diverging i?H should be expected 
[37,38]. Extending the RG approach to the calculation of suitable means for Ru and 
i?L we show that the RG approach can in fact reconcile these findings by establishing 
that the most-probable value of Ru remains rather close to h/e^ even for large Rj^, 
but then diverges with the predicted power- law for i?L > 20h/e^. 

Last, the RG approach is also ideally placed to study the changes of the critical 
properties due to truly long-ranged,* so-called macroscopic inhomogeneities. Close 
to the transition the localization length ^ becomes sufficiently large. Then the long- 
ranged disorder can affect the character of the divergence of ^. We note that in 
previous considerations inhomogeneities were incorporated into the theory through 
a spatial variation of the local resistivity [40^4] , i.e. , "inhomogeneously incoherent . 
Our present approach is able to retain the full quantum coherence. 

The results of this study are presented as follows. The main tool of this work, 
the RG approach to the CC model, is content of section 2. Results for resistance and 
conductance distributions, the level spacing statistics and the associated estimates 
of the critical exponent are reviewed in section 3. The extension of the RG approach 
for the QH insulator and for macroscopic inhomogeneities are presented in section 
4. We conclude in section 5. 



''Note that sometimes the term "long-ranged disorder" is also used for a disorder that has a finite 
correlation radius which is larger than the magnetic length [39]. This is different from the present 
situation. 
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2. The quantum RG approach to the CC model 
2.1. The Chalker-Coddington network model 

The CC network model [12] uses a semi-classical approach to model the integer 
QH transition and is one of the main "tools" for the quantitative study of the QH 
transition [15-17,45-53]. It is based on an extension of the high-field model [54]. In 
order to include the previously mentioned localization-delocalization scenario the 
classical high-field model relies on two basic prerequisites. First, the 2D sample is 
penetrated by a very strong perpendicular magnetic field and second the electrons 
are non-interacting and move in a smoothly varying 2D potential energy landscape 
V{r) illustrated in Fig. 2. As a result the magnetic field B forces an electron onto 



Fig. 2. Left: Illustration of a weakly varying random potential V with equipotential lines at (V) 
and {V) ± A_E which allows to separate the electron motion (black orbit) in a strong magnetic 
field B into cyclotron motion and motion of the guiding center along equipotentials. Right: From 
these equipotentials a,t E = {V), (V) ± AE SPs of the potential can be identified. A single SP acts 
as a scatterer connecting two incoming with two outgoing channels. 

a cyclotron motion with radius Ib much smaller than the potential fluctuations. 
Thus the electron motion can be separated into cyclotron motion and a motion of 
the guiding center along equipotentials of the energy landscape [54] . The cyclotron 
motion leads to the required quantization into discrete Landau levels whereas its 
influence on the electron motion in the potential can be neglected. Besides the 
fundamental assumption about the difference of length scales the model therefore 
contains no further dependence on B. The quantum effects of the electron transport 
through the sample are only determined by the height of the SPs in the potential 
energy landscape. This picture is related to a classical bond percolation problem 
on a square lattice [55] when mapping SPs onto bonds. A bond is connecting only 
when the potential of the corresponding SP equals the potential energy e of the 
electron. From percolation theory [55] follows that an infinite system is conduct- 
ing only at a single e — (V). Consequently, the high-field model describes only a 
single QH transition. It provides a qualitative understanding of the existence of an 
localization-delocalization transition and thus the quantized plateaus in the con- 
ductivity CTxy observed in the integer QH effect [6]. However, because of its purely 
classical assumptions the model is unable to exactly reproduce the critical properties 



<V>+AEA 

<v>\ 

<V>-AE' 




2, 2008 12:32 WSPC/INSTRUCTION FILE review 



6 P. Cain, R. A. Romer 

of the transition found in experiments, i.e., the divergence of the correlation length 
at the transition with an exponent of z/ « 2.3 [56]. Instead, it predicts u = 4/3, the 
value appropriate for classical percolation. 

The CC network model improved the high-field model by introducing quantum 
corrections [12], namely tunneling and interference. Tunneling occurs, in a semi- 
classical view, when electron orbits come so close to each other that the cyclotron 
orbits overlap. From Fig. 2 one can conclude that this happens at the SPs, which 
now act as quantum scatterers described by a unitary scattering matrix S 



(3) 



which connects two incoming with two outgoing channels. Assuming a symmetric 
potential at the SP the scattering rates are given by a pair of a complex transmission 

and a reflection coefficient t and r, respectively. From the required unitarity of S 
it follows that |tp + |rp = 1. Obviously, t and r depend on the potential energy of 
the SP. It was shown [19] that t and r can be parametrized by 

'=(^^)' ''°*' ^'i^y 

where z corresponds to a dimensionless energy difference between SP potential and 
electron energy e. Without restricting the generality {V) = is assumed in the 
following. In case of £ = the value of z then coincides with the dimensionless 
SP height. Similar to bond percolation now a network of SPs can be constructed. 
The SPs are mapped onto nodes and the equipotentials correspond to links. While 
moving along an equipotential an electron accumulates a random phase $ which 
reflects the randomness of the potential. The corresponding phase factor e'* can be 
included in the matrix S [57] or taken into account by additional diagonal matrices 
on each link [19]. 

As in the high-field model this quantum percolation model describes only a 
single QH transition with exactly one extended state in the middle of the Landau 
band at £ = . The critical properties at the transition, especially the value of 
the exponent « 2.4 ± 0.2 [16], agree with experiments [34,56] as well as with 
results of other theoretical approaches [58,59]. For numerical investigations of the 
CC model, one constructs a regular 2D lattice out of the SPs. Then the 2D plane is 
cut into ID slices with the associated scattering matrices transformed into a transfer 
matrix. The conductivity may be calculated by transversing perpendicular to the 
slices along the sample by transfer matrix multiplications [16] according to the 
Landauer-Biittikcr approach [60] . The spatial extension of the 2D plane is limited 
by the computational effort although an additional disorder averaging over many 
samples is not necessary for quasi- ID samples [16]. 

The CC model is a strong-magnetic-field (chiral) limit of a general network 
model, first introduced by Shapiro [61] and later utilized for the study of localization- 
delocalization transitions within different universality classes [62-67]. In addition 
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to the QH transition, the CC model apphes to a much broader class of critical 
phenomena since the correspondence between the CC model and thermodynamic, 
field-theory and Dirac-fermions models [68-76] was demonstrated. 

2.2. Derivation of the basic RG equation for transmission 
amplitudes 

The real-space RG approach [18, 19] can be applied to the CC network analogously 
to the case of 2D bond percolation [77-79]. An RG unit is constructed containing 
several SPs from a CC network. For these SPs the RG transformation has to relate 
their S matrices with the S matrix of the super-SP. The RG unit used here is 
extracted from a CC network on a regular 2D square lattice. The supcr-SP consists 
of five original SPs connected according to Fig. 3. Circles correspond to SPs and 
lines to links in the network. Using this intuitive picture one can identify the loss 
of connectivity in comparison with the original CC network, namely, the four edge 
nodes within a 3 x 3 SP pattern arc fully neglected as are their outer bonds. Thus the 
super-SP has the same number of incoming and outgoing channels as an original SP. 
In analogy to bond percolation the size of the RG unit in terms of lattice spacings 
equals 2 [80]. 




Fig. 3. CC network on a square lattice consisting of nodes (circles) and links (arrows). The RG 
unit combines five nodes (full circles) by neglecting some connectivity (dashed circles). <I>i, . . . , <I>4 
arc the phases acquired by an electron along the loops as indicated by the arrows. ^'i,...,>I'4 
represent wave function amplitudes, and the thin dashed lines illustrate the boundary conditions 
used for the computation of level statistics. 

Between the SPs of the RG unit an electron travels along equipotential lines, and 
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accumulates a certain Aharonov-Bohm phase as in the original network. Different 
phases are uncorrelated, which reflects the randomness of the original potential 

landscape. As shown in section 2.1 each SP is described by an S matrix which 
contributes two equations relating the wave- function amplitudes of incoming Jj , I* 
and outgoing Oi,0* channels. All amplitudes besides the external /i and /| 
can then be expressed by Oi,0* using the phases, e.g. I5 = e'^'^^-'Oi, where $15 
is the phase shift along the link between SPs / and V. The resulting ten modified 
scattering equations form a linear system which has to solved in order to obtain the 
transmission properties of the corresponding super-SP: 



Ax = h 



(5) 
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(6) 



x = iOi, 01 02,0;, 03,0;, 0^,01, 05,01) 



(7) 



and 



b = (ti7i , n/i , 0, 0, 0, 0, r4lt , tJl , 0, 0)' 



(8) 



Note that the amplitudes on the external links coincide with the amplitudes of the 



supcr-SP as h = /* = /'*, O5 



O' and O* = O'*. 



Setting the incoming links of 



the super-SP according to /' = 1, /'* = one can deduce the transmission coefficient 
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t' of the super-SP, since O' = t'l' =t'l = t'. For the transmission coefficient of the 
super-SP this method yields the following expression [19] : 



tit5{r2r3rie^^' - 1) + t2^4e'(*^+^^)(rir3r5e-'^^ - 1) + ^3(^2^56'^^ + ^1^46'^^) 
(rs - r2r4e'*3)(r3 - rirse'*^) + (is - i4t5e'**)(i3 - iii2e**i) 



Here $j corresponds to the sum over the three phases forming a closed loop within 
the RG unit (sec Fig. 3). Equation (9) is the RG transformation, which allows one 
to generate the distribution P{t') of the transmission coefficients of super-SPs using 
the distribution P{t) of the transmission coefficients of the original SPs. Since the 
transmission coefficients of the original SPs depend on the electron energy e, the 
fact that delocalization occurs at e = implies that a certain distribution, Pc{t) — 
with Pc{t^) being symmetric with respect to t"^ = ^ — is the FP distribution of the 
RG transformation (9). 

A systematic improvement of the RG structiire by inclusion of more than five 
SPs into the basic RG unit [57,81,82] leads to similar results. In contrast, using a 
smaller RG unit [38] does not describe the critical properties of the QH transition 
with the same accuracy. The super-SP now consists only of 4 SPs such that it 
resembles the 5SP unit (see Fig. 3) used previously but leaves out the SP in the 
middle of the structure. Again the scattering problem can be formulated as a system 
of now 8 equations which is solved analytically to give 



The result can be verified using Eq. (9) after setting is = and rs = 1, joining 
the phases $1 and $4 and renumbering the indices. As we will show in section 
3.1.3 the value of v obtained for the 4SP RG units is less rehable. Apparently, the 
quality of the RG approach crucially depends on the choice of the RG unit. For 
the construction of a properly chosen RG unit two conflicting aspects have to be 
considered, (i) With the size of the RG unit also the accuracy of the RG approach 
increases since the RG unit can preserve more connectivity of the original network, 
(ii) As a consequence of larger RG units the computational effort for solving the 
scattering problem rises, especially in the case where an analytic solution, as Eq. 
(9), is not attained. Because of these reasons building an RG unit is an optimization 
problem depending mainly on the computational resources available. 

2.3. Description of the RG approach to the LSD 

The nearest-neighbor level-spacing distribution (LSD) is one of the important level 
statistics — other are E2 and A3 statistics — which has been used extensively [83] to 
characterize the universal statistical properties of complex Haniiltonians in the con- 
text of random matrix theory [84]. The necessary eigenenergies are usually obtained 
from the time-independent Schrodinger equation iJ^' = Ek'i' hy diagonalizing the 
Hamiltonian H [85]. After sorting the eigenenergies in ascending order the LSD is 




il/4(/2/;i(.'^'" - l) + /2^i(:"-^'H'■i'■4'•" 



1) 



(10) 



(1 - r2rse'*3)(i _ rir^e^^^) + tit2hUe 



2, 2008 12:32 WSPC/INSTRUCTION FILE review 



10 p. Cain, R. A. Romer 

accumulated from spacings Sk = {Ek+i — -Efc)/A, where Ek+i and are neighbor- 
ing energy levels and A corresponds to the mean level spacing. With the CC model 
based on wave propagation through the sample, H is not accessible directly. In this 
work therefore an alternative approach [86] is used in which the energy levels of 
a 2D CC network can be computed from the energy dependence of the so-called 
network operator U{E). U is constructed similar to Eq. (5). However, when com- 
paring to the calculation of the transmission coefficient t' an essential difference has 
to be taken into account. Energy levels are defined only in a closed system which 
requires to apply appropriate, usually periodic, boundary conditions. The energy 
dependence of U{E) enters through the energy dependence of the ti{E) of the SPs, 
whereas the energy dependence of the phases $j {E) of the links is usually neglected. 
Considering the vector of wave amplitudes on the links of the network, U acts 
similar as a time evolution operator. The eigenenergies can now be obtained from 
the stationary condition 

U{E)^ = (11) 

Nontrivial solutions exist only for discrete energies E^, which coincide with the 

eigenenergies of the system [86]. The eigenvectors correspond to the eigenstates 
on the links. The evaluation of the E^s according to Eq. (11) is numerically very 
expensive. For that reason a simplification was proposed [15]. Instead of solving 
the real eigenvalue problem calculating a spectrum of quasienergies uj is suggested 
following from 

U{E)^i = e''"'^^)*;. (12) 

For fixed energy E the uji are expected to obey the same statistics as the real 
eigenenergies [15]. This approach makes it perfectly suited for large-size numerical 
simulations, e.g. studying 50 x 50 SP networks. 

In order to combine the above algorithm with the RG iteration, in which a 
rather small unit of SPs is considered, some adjustment is necessary. First, one has 
to "close" the RG unit at each RG step in order to discretize the energy levels. From 
the possible variants the closing is chosen as shown in Fig. 3 with dashed lines. For 
a given closed RG unit with a fixed set of tj-values at the nodes, the positions of the 
energy levels are determined by the energy dependences ^j{E) of the four phases 
along the loops. These phases change by '--^ tt within a very narrow energy interval, 
inversely proportional to the sample size. Within this interval the change of the 
transmission coefficients is negligibly small. The closed RG unit in Fig. 3 contains 
10 links and thus it is described by 10 amplitudes. Each link is characterized by 
an individual phase. On the other hand, it is obvious that the energy levels are 
determined only by the phases along the loops. One way to derive U is to combine 
the individual phases into phases $j connected to the four inner loops of the unit. 
The $j are associated with the corresponding "boundary" amplitudes (see Fig. 
3). The original system of 10 equations, which resembles Eq. (5) except for the 
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boundary conditions, can then be transformed to 4 equations by expressing all 
amplitudes in terms of the . The resulting network operator takes the form 



U = 



I (rir2 - ht2t3)e-"^' 

-htirae-''^^ 
\-{t2ri+ht3r2)e-'^^ 



(r4r5 - hUt^)e-'''^^ 



(tir2+i2^3ri)e-'*i 
rir3r4e~'*^ 

-{tit2-hrir2)e-''^^ 

t2r3r5e~'"^^ 

-{hri+t3tAn)e-'-'^^ 



(13) 



Ur2r3e 



-i*3 



r2r-3r5e" 



which can be substituted in Eq. (12). Then the energy levels of the closed RG 
unit including phases ^j{E) = ^j{Ek), are the energies for which one of the four 
eigenvalues of the matrix U is equal to one, which corresponds to the condition 
Lo{Ek) = 0. Thus, the calculation of the energy levels reduces to a diagonalization of 
the 4x4 matrix. It should be emphasized that the reduced size of U in comparison 
with the 10 x 10 matrix resulting from the "straightforward" approach described in 
the beginning directly follows from considering only the relevant energy dependence 
in the four phases of the RG unit. Therefore a larger size of U would lead to 
redundant information for the energy levels. 



2.4. Computing resistances from the transmission coefficients 

We first establish how to connect the experimentally relevant conductances and 

resistances to the transmission coefficients. The dimensionless conductance is simply 
G = The dimensionless longitudinal resistance i?L can be computed from G via 

- ^ - - ^ - 1 = i?2t - 1 (14) 

with the dimensionless 2-terminal resistance i?2t- 

The computation of Ru is less straightforward. We calculate the "resistance" 
R = U/J defined by the potential difference U across the RG unit and the current 

J flowing through the unit. In Fig. 4 this ansatz is illustrated for the 5SP RG imit 
used previously. Assuming that the current enters the RG unit via one incoming 
link (/') only and the other incoming link is inactive (/'* = 0) the resulting power 
transmission coefficient t'^ = 0'^/7'^ of the RG unit can be associated with J. In 
order to determine U one considers the quantities a„ = 0\/I'^ and a,, = /I'"^ as 
chemical potentials measured by weakly coupled voltage probes at these opposite 
links of the RG unit [87]. Thus the voltage drop is given by ?7 = a„ — a^. Because 
of the four-terminal geometry the obtained R contains beside the Hall resistance 
i?H also a contribution from the longitudinal resistance i?L- The separation of the 
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Fig. 4. To determine the Hall resistance Rjj of the RG unit as illustrated here for the 5SP case 
the voltage U and the current J are calculated for positive magnetic field B (left) and the opposite 
direction —B (right). For —B the value of Rjj can be obtained by measuring J either between the 
upper (H-flip) or lower (V-fiip) two links, respectively. 



unwanted part i?L is accomplished by employing the antisymmetry of the HaU 
voltage C/h under the reversal of the direction of the magnetic field B 



Considering the 5SP RG unit again one therefore obtains 



A-B) J-B) 



t' 



12 



(15) 



(16) 



The quantities a„ and a„ ' are calculated by solving the system of equations (5) 
analogously to the determination of t' in section 2.1. One finds 



(ra - r2r4e»*0('^3 - rirse'*^ + {h 



and 



(r3 - r2r4e«*3)(^3 _ nrse**^) + (^3 - Uhe-^^){t^ - tit2e**0 



(17) 



(18) 



Under the reversed field —B the electrons travel along the same equipotentials but in 
the opposite direction. In the corresponding RG unit the links therefore only change 
their direction, as shown in Fig. 4. A full rederivation of the RG equation for —B 
can be omitted if the structure of the RG unit is taken into account. The result with 
identical directions of the current is obtained by flipping the unit vertically (V-flip) 
while the other case corresponds to a horizontal flip (H-flip). The comparison to the 
original RG unit then shows how to map the indices of the SPs and the phases in 
order to adopt the RG equation to — _B, which is demonstrated here for the H-flip 
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result where (ri,ti) ^ {r^,t<^), (r2,<2) ^ (^'4,^4), and $1 ^ $4, 



2 




(19) 



2 




.(20) 



Using Eq. (16) one is now able to determine the distributions of P{Rii) at the QH 
transition iteratively in course of the RG iterations. 

2.5. Iterating the RG structures 

In order to find the FP distribution Pc{t), the RG is started from a certain initial 
distribution of transmission coefficients, Po{t). The distribution is typically dis- 
cretized in at least 1000 bins, such that the bin width is typically 0.001 for the 
interval t G [0, 1]. From Po{t), the ti, i — 1, . . . , 5, are obtained and substituted into 
the RG transformation (9). The phases j = 1, . . . , 4 are chosen randomly from 
the interval <I>j e [0, 27r). In this way at least 10^ super-transmission coefficients t' 
are calculated. In order to decrease statistical fiuctuations the obtained histogram 
Pi{t') is then smoothed using a Savitzky-Golay filter [85]. At the next step the 
procedure is repeated using Pi as an initial distribution. The convergence of the it- 
eration process is assumed when the mean-square deviation J dt[Pn{t) — Pn-i{t)]^ 
of the distribution F„ and its predecessor P„_i deviate by less than 10"'*. 

The actual initial distributions, Po{t), were chosen in such a way that the cor- 
responding conductance distributions, Po(G'), were either uniform or parabolic, or 
identical to the FP distribution found semianalytically [19]. All these distributions 
are symmetric with respect to G = 0.5. One can observe that, regardless of the 
choice of the initial distribution, after 5-10 steps the RG procedure converges to 
the same FP distribution which remains unchanged for another 4-6 RG steps. Small 
deviations from the symmetry with respect to G = 0.5 finally accumulate due to 
numerical instabilities in the RG procedure, so that typically after 15-20 itera- 
tions the distribution becomes unstable and flows toward one of the classical FPs 
P(G) = 5(G) or P(G) = S{G - 1). Therefore the symmetry of the Po(G) with re- 
spect to G = 0.5 is an important requirement in order to converge to the quantum 
FP at all. Note that the FP distribution can be stabilized by forcing P„(G) to be 
symmetric with respect to G = 0.5 in the course of the RG procedure. 

Since the dimensionless SP height Zj and the transmission coefficient U at e = 
are related by Eq. (4), transformation (9) also determines the height of a super-SP 
by the heights of the five constituting SPs. Correspondingly, the distribution P(G) 
determines also a distribution Q{z) of the SP heights via Q{z) = P{G)\dG/dz\ = 
I cosh~^(z/2)P [(e^ -|- 1)^^]- This represents a convenient parameterization of the 
conductance distribution, particularly if small and large t values become important. 
In this case, it is numerically better to perform the RG approach using the Q{z) 
distribution. We typically discretize the distribution Q{z) in at least 6000 bins such 
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that the bin width is 0.01. Since ^; e] — oo, oo[, we have to include lower and upper 
cut-off SP heights such that z e [^;iow, ^up]- Prom Qo{z), we obtain Zi, i = 1, . . . , 5, 

compute the associated t; and substitute them into the RG transformation (9). At 
the next step we repeat the procedure using the new Qi as an initial distribution. 
We check that the values of ziow.up do not influence our results. 

3. Results of the RG procedure at the QH treinsition 
3.1. Fixed point distributions and the critical exponent 
3.1.1. The fixed point distributions Pc{G), Qc{z) 

The distribution Pc{G) of the dimensionless conductance G can be obtained from 
P{t) so that 

Pe(G)^lPe(t) . (21) 

Figure 5 illustrates the RG evolution of P{G) and Q{z). In order to reduce statistical 
fluctuations the FP distribution is averaged over nine results obtained from three 
different Po(G)'s. The FP distribution Pc{G) exhibits a flat minimum around G = 
0.5, and sharp peaks close to G = and G = 1. It is symmetric with respect to 
G « 0.5 with (G) = 0.498 ± 0.004, where the error is the standard error of the 
mean of the obtained FP distribution. Consequently, the FP distribution Qc{z) is 
symmetric with respect to z = 0, which corresponds to the center of the Landau 
band. The shape of Qc{z) is close to Gaussian. We note that the present results agree 
with [57,81,82] where a similar RG treatment of the CC model was carried out. Our 
numerical data have a higher resolution, and show significantly less statistical noise 
because of the faster computation using the analytical solution (9) of the RG [19]. 




Fig. 5. Left; P{G) (thin lines) as function of conductance G at a QH plateau-to-plateau transition. 
Symbols marlf every 20th data point for different initial distributions (■,♦,!), the FP distribution 
(O) e'lid a distribution for RG step n = 16 The vertical dashed line indicates the average of 
the FP distribution. Right; Corresponding plots for the distribution Q{z) of SP heights. 
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3.1.2. The fixed point distributions P{Ru), P{Rl) 

In order to accurately construct the tails of the resistance distributions P(i?L), 
P{Rii), we employ the numerical strategy of section 2.5 and discretize Q{z) with 
ziow = ^20, = 40 for perturbations towards positive z and vice versa for negative 
z. This corresponds to transmission amplitudes t{ziovj) ~ 10~^, t(zup) ~ 1 — 10~^ 
for positive perturbations and vice versa for negative perturbations. Next, we con- 
struct via the RG procedure the sequence of distribution Qn{z), n = 1,2,3,..., 
calciilating at least 10*^ super-transmission coefficients t' and associated z' . We as- 
sume that the iteration process has converged when the mean square deviation 
/ dz [Qn{z) — Qn-i{z)]'^ of the distribution Qn and its predecessor Qn-i deviate by 
less than 10~^. We check that the values of ziow.up do not influence our results and 
that all the previous results at the QH transition as in Refs. [80] are reproduced. 
The full width at half maximum of the FP distribution Qc{z) is about 5 [80]. 

Using Eq. (16) we now obtain besides the FP distribution Qc{z) also the FP 
distribution PdRu)- Like z the value of i?H is not bound to a fixed interval. To 
perform a discretization requires to set lower and upper bounds appropriately. The 
resulting histogram is limited to [—100,100] containing 40000 bins. Since Pc{t) is 
already known from section 3.1.1 one can speed up the determination by using Qc{z) 
as initial Qo{z) and therefore obtain Pc{Rh) already after the first iteration. 

In Figure 6 the resulting distributions Pc(-Rl) and Pc{Rh) have been plotted. 
Both distributions are strongly non-Gaussian with long tails. Pc{Rh) can be fitted 
by a log-normal distribution [37]. We have computed Pc{Rb_) using both V- and 
H-flip structures. Their graphs in Fig. 6 show a significant difference in the shape of 
Pc{Rh)- The H-flip distributions are characterized by a very sharp, nearly symmetric 
peak at Rh = 1 which coincides with the value of Rh at the first Landau level. For 
V-flip an asymmetric distribution is found again with a strong peak at Ru = 1- 
Surprisingly, there exists also a kink at Ru = 0.5. Such a peak would normally 
correspond to a second Landau level, which is however not part of the CC model. 
The origin of the "wrong" kink must stem from a difference between the calculation 
of Ru according to H-fiip and V-flip method, respectively. Both methods differ only 
in the {—B) part of Ru- As illustrated in Fig. 4, the H-flip method measures J 
between the same SPs / and V as for {+B). In contrast, the J for V-flip is obtained 
between SPs // and IV. While one can expect the same result for H-flip and V-flip 
when considering only the (—B) contribution separately, the flnal V-flip result for 
Ru according to Eq. (16) thus corresponds to a measurement in a different sample 
with different random potential. We therefore only use the H-flip results in the 
following. 

3.1.3. Critical exponent 

The language of the SP heights provides a natural way to extract the critical 

exponent v. Suppose that the RG procedure starts with an initial distribution 
Qo{z) = Qc{z — zq) that is shifted from the critical distribution Qc{z) by a small 
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Fig. 6. Left: Distribution of the longitudinal resistance Kl- The □ symbols indicate the FP 
distribution, the solid line is a fit to a log-normal distribution for positive Hl- Only every 5th data 
point is shown. Right: Distribution of the Hall resistance Rh at the FP as obtained from the V-flip 
(dashed line) and the H-flip method (solid line). Note the kink at Rh = 1 for the V-flip structure. 



zo oc e. The meaning of zq is an additional electron energy measured from the center 
of the Landau band. The fact that the QH transition is infinitely sharp at = im- 

pHes that for any zq =/= 0, the RG procedure drives the initial distribution Q{z — zq) 
away from the FP. Since zq <C 1, the first RG step would yield Qc{z — tzq) with 
some number r independent of zq- At the nth step the center of the distribution 
will be shifted by Zmax,n = t'^zq, while the sample size will be magnified by 2". 
After a certain number of steps, say Hl, the shift will grow to 

•^max, ni. — T ^ Zq 1. (^^) 

where a typical SP is no longer transmittablc. Then the localization length S, can be 
identified with the system size 2"^ a where a is the lattice constant of the original 
RG unit. Using this relation one can rewrite r"^ in Eq. (22) as a power of (C/a) 

(e/a) 2) _ ;L (23) 

from which follows that ^ diverges as 

with u = In 2/ In T. When the RG procedure is carried out numerically, one should 
check that zq is small enough so that Zmax,n oc zq for large enough n. Consequently, 
the working formula for the critical exponent can be presented as 



In 



which should be independent of n for large n. 

Thus, starting with an initial shift of Qc{z) by a value zq results in the further 
drift of the maximum position, Zmax,n, away from z = after each RG step. As 
expected, ZinaK,n depends linearly on zq as shown in Fig. 7 (inset) for different 
n from 1 to 8. The critical exponent is calculated from the slope according to Eq. 
(25). Its value converges with n to 2.39±0.01. The error corresponds to a confidence 
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interval of 95% as obtained from the fit to a linear behavior. In Fig. 7, we also show 
results obtained for the 4SP structure discussed in section 2.2. Obviously, 1/43? tends 
towards a different value. 



2 4 5 6 7 8 




4 16 32 64 128 256 



2" 

Fig. 7. Critical exponent u obtained by the QH-RG approach as function of effective Uncar 
system size L = 2" for RG step n. Solid and open symbols correspond to the 5SP and 4SP RG 
unit, respectively. The error bars correspond to the error of linear fits to the data. The dashed line 
shows u = 2.39. Inset: 1/ (here for 5SP) is determined by the dependence of the maximum 2niax,n 
of Qn{z) on a small initial shift zq. Symbols indicate the eight RG steps in accordance with the 
main plot. Dashed lines indicate the linear fits. 



3.1.4. Comparison with other works 

By dividing the CC network into units, the RG approach completely disregards at 
each RG step the interference of the wave-function amplitudes between different 
units. For this reason it is not clear to what extent this approach captures the 

main features and reproduces the quantitative predictions at the QH transition. 
Therefore, a comparison of the RG results with the results of direct simulations 
of the CC model is crucial. These direct simulations are usually carried out by 
employing either the quasi-lD version [12,16,88] or the 2D version [89-92] of the 
transfer-matrix method. For the critical exponent the values v = 2.5 ± 0.5 [12] and 
later i/ = 2.4 ± 0.2 [16] were obtained. Note that the result of the previous section 
is in excellent agreement with these values, and is also consistent with the most 
precise jy = 2.35 ± 0.03 [14]. This already indicates the remarkable accuracy of the 
RG approach using the 5SP structure. 
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In [93] and [94] Pc{G) was studied by methods which are not based on the 
CC model. Both works reported a broad distribution Pc{G). In [90-92] the critical 

distribution Pc{G) of the conductance was studied. Pc{G) was found to be broad, 
which is in accordance with Fig. 5. One can compare the moments of Pc{G) to those 
calculated in [90,92]. They agree with the present RG calculations up to the sixth 
moment. We emphasize that the moments from [92] can hardly be distinguished 
from the moments of a uniform distribution. This reflects the fact that Pc{G) is 
practically flat except for the peaks close to G = and G = 1. 

3.2. Energy-level statistics at the QH transition 

3.2.1. Choosing an appropriate energy dependence of the phases 

The crucial ingredient for obtaining the LSD is the choice of the energy dependence 
^j{E). If each loop in Fig. 3 is viewed as a closed equipotential as it is the case for 
the first step of the RG procedure [12], then ^j{E) is a true magnetic phase, which 
changes linearly with energy with a slope governed by the actual potential profile 
which in turn determines the drift velocity. Thus we make the ansatz 

= $0 J +2^-, (26) 

where a random part <I>o.j is uniformly distributed within [0,27r), and 27r/,Sj is a 
random slope. Here the coefficient Sj acts as an initial level spacing connected to 
the loop j of the RG unit by defining a periodicity of the corresponding phase. 
Strictly speaking, the dependence (26) applies only for the first RG step. At each 
step n > 1, ^j{E) is a complicated function of E which carries information about 
all energy scales at previous steps. However, in the spirit of the RG approach, one 
can assume that ^j{E) can still be linearized within a relevant energy interval. 
The conventional RG approach suggests that different scales in real space can be 
decoupled. Linearization of Eq. (26) implies a similar decoupling in energy space. 
In the case of phases, a "justification" of such a decoupling is that at each RG step, 
the relevant energy scale, that is the mean level spacing, reduces by a factor of four. 

With ^j{E) given by Eq. (26), the statistics of energy levels determined by 
the matrix equation (12) is obtained by averaging over the random initial phases 
$o,j and values tt chosen randomly according to a distribution P{t). For every 
realization the levels Ek are computed from the solutions uj{Ek) = of Eq. (12) 
as illustrated in Ref. [83]. The energy interval is scanned in discrete energy steps 
with AE = min{sj}/250 which is adapted to each random realization of the $oj 
and takes the periodicity in Eq. (26) and its influcince on the behavior of u>{E) into 
account. In particular, each realization yields three level spacings which are then 
used to construct a smooth LSD. Thus the situation is comparable with estimating 
the true RMT ensemble distribution functions from small, say, 2x2 matrices only 
[84,95]. The outline of the RG procedure for the LSD is as follows. The slopes 
Sj in Eq. (26) determine the level spacings at the first step. They are randomly 
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distributed with a distribution function Po{s). Subsequent averaging over many 
realizations yields the LSD, Pi(s), at the second step. Then the key element of the 
RG procedure, as applied to the level statistics, is using Pi(s) as a distribution of 
slopes in Eq. (26). This leads to the next-step LSD and so on. 

The approach of this work relies on the "real" eigenenergies of the RG unit. 
The simpler computation of the spectrum of quasienergies adopted in large-scale 
simulations within the CC model [15,17] cannot be applied since the energy depen- 
dence of phases $ j in the elements of the matrix is neglected and only the random 
contributions, ^o,j, are kept. Nevertheless it is instructive to compare the two pro- 
cedures. In Rcf. [83] we have calculated the dependence of the four quasienergies 
Uk on the energy E with ti chosen from the critical distribution Pc{t). The energy 
dependence of the phases i>j was chosen from LSD of the Gaussian unitary ensem- 
ble (GUE) [84,95] according to Eq. (26). We find that the dependences u){E) range 
from remarkably linear and almost parallel to strongly nonlinear. 

3.2.2. The critical LSD at the QH transition 

Let us first consider the shape of the LSD at the QH transition. As starting dis- 
tribution Po{s) of the RG iteration we choose the LSD of GUE because previous 
simulations [13, 15] indicate that the critical LSD is close to GUE. According to 
Po(s), each Sj is selected randomly and <I>j, j = 1, ... ,4 is set as in Eq. (26). For 
the transmission coefficients of the SP the FP distribution (t) is used. As known 
from section 3.1.1, Pn{t) drifts away from the FP within several further iterations 
due to unavoidable numerical inaccuracies. In order to stabilize the calculation, 
the FP distribution Pc{t) is therefore used in every RG step instead of Pn{t)- This 
trick does not alter the results but speeds up the convergence of the RG for Pc{s) 
considerably. By finding solutions w{Ek) = of Eq. (12) the new LSD Pi(s') is 
constructed from the "imfolded" energy level spacings s'f, = (Ek+i — Ek)/A, where 
k = 1,2,3, and the mean spacing A = {E4 — Ei)/3. Due to the "unfolding" [96] 
with A, the average spacing is set to one for each sample and in each RG-iteration 
step spacing data of 2 x 10^ super-SPs are taken into account. The resulting LSD 
is discretized in bins with largest width 0.01. In the following iteration step the 
procedure is repeated using Pi (s) as initial distribution. Convergence of the itera- 
tion process is assumed when the mean-square deviation of P„(s) deviates by less 
than 10~* from its predecessor P„_i(s). The above approach now enables one to 
determine the critical LSD Pc(s). The RG iteration converges rather quickly after 
only 2 — 3 RG steps. The resulting Pc(s) is shown in Fig. 8 together with the LSD 
for GUE. 

Although Pc(s) exhibits the expected features, namely, level repulsion for small 
s and a long tail at large s, the overall shape of Pc(s) differs noticeably from GUE. 
In the previous large-size lattice simulations [13, 15] the obtained critical LSD was 
much closer to GUE. This fact, however, does not necessarily imply a lower accu- 
racy of the RG approach. Indeed, as it was demonstrated recently, the critical LSD 
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Fig. 8. Left: FP distributions Pc(s) obtained from the spectrum of uji (E = 0) and from the RG ap- 
proach using the real eigenenergies Ek in comparison to the LSD Pgue(s) = 32s'^ exp {—^s'^)/tt'^ 
for GUE. As in all other graphs P{s) is shown in units of the mean level spacing A. Right: Pcis) 
for small s in agreement with the predicted behavior. Due to the log-log plot errors are shown 
in the upper direction only. 

- although being system-size independent — nevertheless depends on the geometry 
of the samples [97] and on the specific choice of boundary conditions [98,99]. Sen- 
sitivity to the boundary conditions does not affect the asymptotics of the critical 
distribution, but rather manifests itself in the shape of the "body" of the LSD. One 
should note that the boundary conditions which have been imposed to calctilate the 
energy levels (dashed lines in Fig. 3) are non-periodic in contrast to usual large-size 
lattice simulations [15]. 

As mentioned above there is also the possibility to assess the critical LSD via 
the distribution of quasienergies. In Fig. 8 the result of this procedure is shown. It 
appears that the resulting distribution is almost identical to Pc{s). This observation 
is highly non-trivial, since there is no simple relation between the energies and 
quasienergies [83]. Moreover, choosing another functional form for ^{E) instead of 
the linear i?-dcpcndcnce (26), the RG procedure yields an LSD which is markedly 
difl:erent (within the "body") from Pc{s) [83]. Using the quasienergies instead of 
real energies, as in [15], and linearizing the energy dependence of phases [as in Eq. 
(26)] is not rigorous. However, the coincidence of the results for the two procedures 
supports the validity of the approach. 

3.2.3. Smalls and large-s behavior 

The general shape of the critical LSD is not universal. However, the small-s behavior 
of Pc{s) must be the same as for GUE, namely Pc{s) oc . This is because delo- 
calization at the QH transition implies level repulsion [100,101]. Earlier large-scale 
simulations of the critical LSD [13,15,17,102-109] satisfy this general requirement. 
The same holds also for the result of this work, as can be seen in Fig. 8. The given 
error bars of the numerical data are standard deviations computed from a statistical 
average of 100 FP distributions each obtained for different random sets of tj's and 
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$j's within the RG unit. In general, within the RG approach, the s^-asyniptotics 
of P{s) is most natural. This is because the levels are found from diagonalization 
of the 4x4 unitary matrix (13) with absolute values of elements widely distributed 
between and 1. 

The right form of the large-s tail of P(.s) is Poissonian with Pc{s) oc exp(— 6s) [13, 
101,103,110]. It appears that the accuracy of the present RG approach is insufficient 
to discern this non-trivial feature of the critical LSD, since we have a good accuracy 
only for s/A < 2.5 [83]. 

3.2.4. Finite-size scaling at the QH transition 

In order to extract u from the LSD the one-parameter-scaling hypothesis [111] is 
employed. The approach describes the rescaling of a quantity a{N-, {zi}) - depend- 
ing on (external) system parameters {zi} and the system size — onto a single 
curve by using a scaling function / 



We now use the scaling assumption to extrapolate f to N ^ oo from the finite-size 
results of the compntations. The knowledge about / and ^oo then allows to derive 
the value of ;/ similarly to (24). In order to define a suitable control parameter zq 
in the transition region, again the natural parameterization (4) of the transmission 
coefficients is used. Because the QH transition occurs exactly at i;o = any initial 
distribution Q(z — zq) with zq =^ will evolve during the RG procedure toward an 
insulator, either with complete transmission of the network nodes (for zq > 0) or 
with complete reflection of the nodes (for zq < 0). 

In principle, one is free to choose for the finite-size scaling (FSS) analysis any 
characteristic quantity a{N; za) constructed from the LSD which has a systematic 
dependence on system size N for zq ^ while being constant at the transition 
zo = 0. Out of the large number of possible choices [13,101,110,112-114] a restriction 
is made to quantities that are defined mainly by the small-,s behavior which is 
accurately described by the RG approach. The quantities are obtained by integration 
of the LSD and have already been successfully used in [112, 113, 115], namely 



with I{s) = Jq P{s')ds'. The integration limit for the first 2 quantities is chosen as 
So = 1.4 which approximates the common crossing point [112] of all LSD curves as 
can be seen in Fig. 9. Thus P{so) is independent of the distance \z — Zc\ to the critical 
point and the system size N. Note that N is directly related to the RG step n by 
N = 2". The double integration in ai is numerically advantageous since fluctuations 
in P{s) are smoothed. One can now apply the finite-size-scaling approach from Eq. 
(27) for ap^i^s{N,zo)- Since Q;p,i,s(-^5 ^o) is analytical for finite N, one can expand 




(27) 




(28) 
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Fig. 9. RG of the LSD used for the computation of i/. The dotted hnes corresponds to the first 9 
RG iterations with an initial distribution Pq shifted to the metalhc regime (zq = 0-1) while the full 
thin lines represent results for a shift toward localization (zo = —0.1). Within the RG procedure 
the LSD moves away from the FP as indicated by the arrows. At s/A R! 1.4 the curves cross at 
the same point which is used when deriving a scaling quantity from the LSD. 

the scaling function / at the critical point. The first-order approximation yields 

a{N, zq) ~ q;(A'', Zc) + a\zo — Zc IN^/" (29) 

where a is a coefficient. Better results are obtained using a higher-order expansion 
proposed by Slevin and Ohtsuki [116], but contributions from an irrelevant scaling 

variable can be neglected since the transition point zq = is known [83]. 

In order to obtain the functional form of / the fitting parameters, including v, 
are evaluated by a nonlinear least-square (x^) optimization. In Fig. 10 the resulting 
fits for, e.g., as at the transition arc shown. The fits are chosen such that the total 
number of parameters is kept small and the fit agrees well with the numerical data. 
The corresponding scaling curve is displayed in the right panel of Fig. 10. In the 
plots the two branches for complete reflection (^o < 0) and complete transmission 
{zq > 0) can be distinguished clearly. In order to estimate the error of the fitting 
procedure the over 100 results for v obtained by different orders of the expansion, 
system sizes A'', and z ranges around the transition are compared [83,117]. The value 
oiv = 2.37 ± 0.02 is calculated as average over many individual fits for apj. It is in 
excellent agreement with 2.39 ±0.01 calculated in section 3.1.3. The value obtained 
from the parameter-free scaling quantity as = s^P{s)ds [113] is = 2.33 ±0.05, 
but scaling is less reliable due to the insufficient statistics in the large-s tail. 

Finally the influence of the initial conditions on the result of the LSD and 
the one-parameter scaling has been studied in Ref. [83]. The corresponding LSD 
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Fig. 10. Behavior of ag at the QH transition. Left: Data for RG iterations n = 1, ... ,9 corre- 
sponding to effective system sizes JV = 2" = 2, . . . , 512. Full lines indicate the functional depen- 
dence according to FSS using the minimization with Oi, = 2 and Oz = 3. Because of large 
deviations data for AT = 2 were excluded from the fitting. Right: Resulting FSS curves. Different 
symbols correspond to different effective system sizes AT = 2". 



agrees with GUE less convincingly than the LSD computed using the true Pc(G). 
Estimates of the critical exponent yield also less accurate values vi = 2.43 ± 0.02 
and = 2.46 ±0.03 which are nevertheless still reasonably close tov = 2.37 ±0.02. 
Thus our choice of the initial distribution Pc{G) is indeed valid. 



4. Extensions of the RG approach 
4.1. The QH-to-insulator transition 

4.1.1. Evolution of the resistance distributions P{Ri^) and P{RyC} 

In order to model the transition into the insulating regime, we shift the initial 
distribution Qo{z) Qo{z — zo) by a small zq. For zo < and zo > 0, the RG 
flow will then drive the distributions into the insulating, G ^ 0, and plateau, 
G — > 1, regimes, respectively. In Fig. 11, we show the resulting distributions after 
many RG steps. For P{Ri^), we see in Fig. 11 that the flow towards G ^ 1 results 
in a decrease of large J?l events, whereas conversely, the regime G — > leads to 
an increase in large i?L values and a decrease of the maximum value in P{Ri^). For 
P{Rh), Fig. 11 shows that the plateau regime G — > 1 gives a highly singular peak at 
the dimensionless quantized Hall value 1, corresponding to a perfect Hall plateau. 
On the other hand, the insulating regime G ^ shows an increase of weight in 
the tails of P(i?H) and the eventual obliteration of any central peak. Very large 
absolute values of i?L and i?H in the tails are obviously the consequence of a very 
small denominator t'^ in Eqs. (14) and (16) meaning that the current is almost 
totally reflected by the RG unit. This scenario is possible since P„(t^) also exhibits 
strong fluctuations and spreads over the whole range between and 1. 
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Fig. 11. Left: Distribution of the longitudinal resistance Rj^. The □ symbols indicate the FP 
distribution, the solid line is a fit to a log-normal distribution and the dot-dashed and dashed lines 
show the distribution after n = 15 RG iterations into the conductance regimes G{zo < 0) ^ 
and G{zo > 0) — > 1. Only every 5th data point is shown for P(Rl). Right: Distribution of the 
Hall resistance Rh- The FP distribution has been shaded to zero. For the G{zo < 0) ^ and 
G{zo > 0) — > 1 regimes, we have indicated the distributions after n = 16 RG iterations by bold 
dashed and dot-dashed lines. 



4.1.2. Evolution of the average resistances 

In order to extract an averaged and Ru from these non-standard distributions, 
we should now select an appropriate mean (•) that characterizes and captures the es- 
sential physics and allows comparison with the experimental data. The precise oper- 
ational definition is also important as it corresponds to different possible experimen- 
tal setups. Therefore, we consider several means: (i) arithmetic (i?)ari = 
(ii) geometric/typical (i?)typ = cxp^^ lni?i/A'^, (iii) median (central value) (i?)mod, 
where N denotes the number of samples (> 10^) in each case. The median and the 
typical mean (and their variances [38]) are less sensitive to extreme values than other 
means (such as, e.g. root-mean-square and harmonic mean) and this makes them a 
better measure for highly skewed and long-tailed distributions such as P{Rh) and 
P{Rh) in the insulating regime. In the plateau regimes, the distributions are less 
skewed, particularly for Ru and the difference in the means becomes less important. 

We arc left with determining in which operational order to apply the averaging 
procedure. For R^, as measured via Ohm's law (14) as a ratio, it is obvious that the 
appropriate average should be i?L = — 1 (and not (^) — 1). For Ru, the situation 
is less straightforward due to the definition of J7h in (16). A simple average is {Uji), 
i.e. using the appropriate P(i?H)- Similar to the experimental procedure, we can 
also estimate Uu via (a„ — a^) for each B field direction separately. In Ref. [38], it 
has been suggested that a more appropriate average {Uu)* can be constructed from 
{o-u) — {(iv)- This later procedure corresponds to measuring the voltage drop between 
positions Uu and Uy in Fig. 3 by separately measuring the individual voltages with 
respect to ground and then recombining them. For _Rh , this yields 

1 yxC-u )typ - [av )typ J - \^\au )typ - [av )typ J 

{Ru):y, = 2 . (30) 



I 



I 
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o most-probable value^of P(Rh)' 




- o most-probable value of P{Rh) 
arithmetic <RH>art 
typical <Rh>*i,p 
median of PiR^) xkx 




Fig. 12. Left: Dependence of averaged Rn on averaged i?L for various means obtained from 10* 
samples. The dashed line describes the divergence Ru oc R^ of the typical mean (i?H)typ> the 
median (i?H)med a^iid the most probable value of P{Rji). The variance of the eirithmetic mean 
diverges as Rj^ — » oo and (RH)ari is no longer a useful characteristic of P(ilH)- The horizontal 
dotted line indicates Rh = 1- Right: Plot of the same data as in Fig. 12 but on a linear scale and 
for experimentally accessible resistance values. Here the almost quantized behavior for the most 
probable value of P(Rji) becomes even more pronounced. The horizontal dotted line indicates 
Rh = 1. 



In Fig. 12 we show the resulting dependence /JhC-Rl) when these averaging defini- 
tions are being used. In the plateau regime i?L < 1, the nearly constant behavior 
of all i?H averages is as expected. For i?L oo, we get a divergent i?H when using 
the median and typical means as suggested by Refs. [37, 38] to compute both Ri, 
and i?H- This divergence can be captured by a power-law i?H with k « 0.26. 

The arithmetic mean for large Ri, ^ 1 quickly becomes instable and no useful 
information can be inferred. 

Reducing the information to the experimentally more relevant resistance regime 
of a few times h/e"^, we replot the i?L and i?H data in Fig. 12 (right panel) on 
a linear scale. In addition to the three means above, we also show the behavior 
of the most-probable value Rh at which P{Rh) has a maximum. This estimate 
Rh{Rl) appears relevant in the experimental setup where 10^ different samples 
cannot be easily measured and the full distribution functions cannot be constructed 
in similar detail. Most importantly, for the range of i?L values shown in Fig. 12, the 
value of Rh deviates only slightly from its quantized value 1 at the transition and 
Ri^ + Rh « R2t [36]. Therefore, the experimental estimate of Rh appears to support 
the notion of the quantized Hall insulator. Indeed, the deviations from 1 are less 
than 10% until i?L ~ 25. However, going back to the left panel of Fig. 12, we see 
that in the strongly insulating regime i?L oo also Rh diverges with a power-law 
that is well-described by {R)lyp. We emphasize that fluctuations in _Rh for large 
i?L ^ 10^ are due to numerical inaccuracies in P{Rh) and decrease upon further 
increasing the number of samples. 

The results for i?L and Rh in the localized regime can be very well described 
by an exponential scaling function with finite-size correction [37] i?L,H(2",2) oc 



2, 2008 12:32 WSPC/INSTRUCTION FILE review 



26 P. Cain, R. A. Romer 

2'''" exp[2"^L ^(z)]. Plotting 5l,h as a function of small perturbation zq, we find 
^l,h(2o) ^0 '^^'^ with i/l,h ~ 2.35 as shown in Fig. 13. Thus we recover the univer- 
sal divergence of the localization length ^ even when using the RG in the insulating 
regime [14] . An equally reliable estimate of the irrelevant exponent 7 appears not 




Fig. 13. Left: Dependence of the localization lengths and on the initial shift zq using the 
ansatz (8) of [37]. We obtain a reasonable estimate of the critical exponent u. Right; The semi 
circle law for different means as in Fig. 12. The most probable values (o) show perfect semi-circle 
(solid line) behavior. The thin dotted line denotes ctl = 1/2. 

possible for our data. As usual, the 4-terminal resistances can be converted into 

the respective conductances via ctl.h = Rlm/ {Rl + -^h)- ^'^^ these conductances, 
one expects the semi-circle law + {an — 1/2)^ = (1/2)^ to hold [52]. In Fig. 13, 
we show that the most-probable values capture the overall shape and symmetry 
properties best [25,36], the other averages show pronounced deviations. We also 
note that the relation Ri^(zq) = l/Ri^(—zo) is obeyed by all means [118]. This is a 
consequence of the reflection symmetry of P{G) [80]. 

4.2. Macroscopic inhomogeneities in the RG approach 

4.2.1. Implementing power-law correlations in the RG approach 

A natural way to incorporate a quenched disorder into the CC model is to ascribe 
a certain random shift, zq, to each SP height, and to assume that the shifts at 
different SP positions, r and r', are correlated as 

{zQ{r)zQ{r'))^\r-r'\-'\ (31) 

with a > 0. After this, the conventional transfer-matrix methods of [88, 89] could 
be employed for numerically precise determination of (G), the distribution Pc{G), 
its moments, and most importantly, the critical exponent, v. However, the transfer- 
matrix approach for a 2D sample is usually limited to fairly small cross sections (e.g., 
up to 128 in [91]) due to the numerical complexity of the calculation. Therefore, the 
spatial decay of the power- law correlation by, say, more than an order in magnitude 
is hard to investigate for small a. 
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On the other hand, the RG approach is perfectly suited to study the role of 
the quenched disorder. First, after each step of the RG procedure, the effective 

system size doubles. At the same time, the magnitude of the smooth potential, 
corresponding to the spatial scale r, falls off with r as r~"/^. As a result, the 
modification of the RG procedure due to the presence of the quenched disorder 
reduces to a scaling of the disorder magnitude by a constant factor 2^"/^ at each RG 
step. Second, the RG approach operates with the conductance distribution Pn{G) 
which carries information about all the realizations of the quenched disorder within 
a sample of size 2". This is in contrast to the transfer-matrix approach [88,89], 
within which a small increase of the system size requires the averaging over the 
quenched disorder realizations to be conducted again. 

The above consideration suggests the following algorithm. For the homogeneous 
case all SPs constituting the new supcr-SP are assumed to be identical, which means 
that the distribution of heights, Q,, (z), for all of them is the same. For the correlated 
case these SPs are no longer identical, but rather their heights are randomly shifted 
by the long-ranged potential. In order to incorporate this potential into the RG 
scheme, Qn{zi) should be replaced by Qn{zi — A^"-*) for each SP, i, where aI"-* is 
the random shift. Then the power-law correlation of the quenched disorder enters 
into the RG procedure through the distribution of a\"'\ That is, for each n the 
distribution is Gaussian with the width /3(2")~"/^. For large enough n the critical 
behavior should not depend on the magnitude of the correlation strength /?, but on 
the power, a, only. 

4.2.2. A new critical behavior 

We find that for all values of a > in correlator (31) the FP distribution is identical 
to the uncorrelated case within the accuracy of the calculation. In particular, (G) — 

0.498 is unchanged. However, the convergence to the FP is numerically less stable 
than for uncorrelated disorder due to the correlation- induced broadening of Qn{z) 
during each iteration step. In order to compute the critical exponent v{a) the RG 
procedure is started from Qo{z — zq), as in the uncorrelated case, but in addition 
the random shifts are incorporated. As already explained these shifts are a result 
of the quenched disorder in generating the distribution of z at each RG step. The 
outcome is shown in Fig. 14. It illustrates that for increasing long-ranged character 
of the correlation (decreasing a) the convergence to a limiting value of v slows down 
drastically. Even after eight RG steps (i.e., a magnification of the system size by a 
factor of 256), the value of u with long-ranged correlations still differs appreciably 
from V — 2.39 obtained for the uncorrelated case. The RG procedure becomes 
unstable after eight iterations, i.e., Zmax.g can no longer be obtained reliably from 
Qg{z). Unfortunately, for small a the convergence is too slow to yield the limiting 
value of V after eight steps only. For this reason, the above method is, strictly 
speaking, unable to unambiguously answer the question whether sufficiently long- 
ranged correlations result in an a-dependent critical exponent v{ol), or whether the 
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value of v eventually approaches the uncorrelated value of 2.39. Nevertheless, the 
results in Fig. 14 indicate that the effective critical exponent exhibits a sensitivity 
to the long-ranged correlations even after a large magnification by 256 x 256. Since 
this factor coincides with the change in scale from the microscopic magnetic length 
to the realistic samples with finite sizes of several ^m, macroscopic inhomogeneities 
are able to affect the results of scaling studies. Note further that as shown in the 



Iogio2 




2" 12 3 4 
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Fig. 14. Left: Critical exponent 1/ obtained by the QH-RG approach as a function of RG scale 2" 
for (3 = 3 and different correlation exponents a. The dashed line indicates f = 2.39, which is the 
value obtained for uncorrelated disorder. For clarity, the errors are shown only for a = 0.2 and 4. 
Inset: 1/ vs 2°"/-^ shows deviations from scaling. Right: Dependence of the critical exponent u on 
correlation exponent a for different (3 = 1,2,3, and 4 as obtained after eight QH-RG iterations. 
The dotted line indicates Oc = 0.84. The dot-dashed line u = 2/a for a < ac follows from the 
extended Harris criterion [119] for classical percolation. 

inset of Fig. 14 one can observe scaling of u values when plotted as function of a 
renormalized system size 2""/^ only for large correlation strength /?. In this case 
the disorder broadening during the first RG steps becomes comparable to the FP 
distribution. Therefore in the QH-RG the correlated disorder dominates over the 
initial FP distribution. One should emphasize that ^{a) obtained after eight RG 
steps always exceeds the uncorrelated value. 

Figure 14 shows the values of p obtained after the eighth RG step as a function 
of the correlation exponent a for different dimensionless strengths (3 of the quenched 
disorder. It is seen that in the domain of a, where the values of v differ noticeably 
from 1/ = 2.39, their dependence on (3 is strong. According to the extended Harris 
criterion [119,120], i/(a) is expected to exhibit a cusp at the /3-independent value 
ac = 2/2.39 w 0.84. From the results in Fig. 14, two basic observations can be 
made. For a small enough magnitude of the long-range disorder, one observes a 
smooth enhancement of iy{a) with decreasing a without a cusp. Although such a 
behavior might be caused by the relatively small number of RG steps, the data 
can be relevant for realistic samples which have a finite size and a finite phase- 
breaking length governed by temperature. At the largest /3 = 4, there is still no 
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clear cusp but the a-dependence for a < ac is closer to the v = 2/a predicted by 
the extended Harris criterion. Unfortunately, the numerics becomes progressively 
unstable, forbidding to go to even larger (3. 

4.2.3. Intrinsic short-range disorder in quantum percolation 

Here one should note that there is a crucial distinction between the classical case and 
the quantum regime of the electron motion considered in the present RG approach. 

In this study, the correlation of the heights of the SPs has been incorporated into 
the RG scheme. At the same time it was assumed that the Aharonov-Bohm phases 
acquired by an electron upon traversing the neighboring loops are uncorrelated. This 
asstunption implies that, in addition to the long-ranged potential, a certain short- 
ranged disorder causing a spread in the perimeters of neighboring loops of the order 
of the magnetic length is present in the sample. The consequence of this short-range 
disorder is the sensitivity of the results to the value of /3 which parametrizes the 
magnitude of the correlated potential. The presence of this short-range disorder 
affecting exclusively the Aharonov-Bohm phases significantly complicates the ob- 
servation of the cusp in the i/(a) dependence at a w 0.84, as is expected from the 
extended Harris criterion. A more detailed investigation of this effect has been given 
in Ref. [80]. We emphasize that an unambiguous demonstration of the validity of 
the extended Harris criterion has recently been given by the Liouville approach to 
the random Landau matrix model [121]. 

5. Conclusions 

Real-space RG approaches to problems of quantum statistical physics have been 
applied with much success in many context. The RG necessarily uses only a lo- 
cal approximation to the full connectivity of the network. Thus it is important to 
choose the RG structures such that the overall symmetries of the network are re- 
tained. In the present work, we have shown various examples where such a strategy 
is permissible in the context of the integer QH effect. Particularly for the deter- 
mination of critical distributions Pc{G), Pc{Rh), Pc{Rii) and the critical exponent 
v = 2.39 ± 0.01, this allowed us to study of the critical properties of the QH tran- 
sition with surprisingly high accuracy. It seems possible to further enhance the 
accuracy of the RG by increasing the size of the RG unit [82]. However, in order 
to fully benefit from the increased size, it is essential that an analytic solution as 
in Eq. (9) is obtained. Otherwise, numerical round-off and stability effects will lead 
to an overall decreased accuracy. We emphasize also that in the case of larger RG 
units the size-rescaling factor needs to be calculated with care. 

The robustness of the 5SP RG when computing universal properties is more- 
over demonstrated by an FSS analysis of the LSD around the QH transition. The 
exponent u = 2.37 ± 0.02 of the localization length obtained by a nonlinear 
minimization is in excellent agreement with the value reported in the literature. It 
is surprisingly good when keeping in mind that it was derived just from the four 
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loops of the RG unit which therefore seem to capture the essential physics of the 
QH transition. The success of the RG approach to the LSD can be attributed to 

the description of the transmission amplitudes t of the SPs by a full distribution 
function P{t) while for network models usually a fixed value t{E) is assigned to all 
SPs for simplicity. Also, the phases are associated with full loops in the network and 
not with single SP-SP links. Due to these reasons the RG iteration is always gov- 
erned by the quantum critical point of the QH transition. Even when starting the 
iteration with a distribution Po{G) totally different from Pc{G), but still symmetric 
with respect to G = 0.5, one approaches the same results after a few iterations. 

Having established that the RG approach works at the QH transition, wc then 
set out to extend the method to distributions away from criticality. We have shown 
in a quantum coherent calculation that the insulating Hall behavior of i?H(^?L — * oo) 
is dominated by the power- law divergence i?H oc R^. However, in the experimentally 
attained range i?,L < 10 [25,26,36], the deviation from a quantized i?,H = 1 is very 
small and the onset of the divergence for i?L ^ 10 is yet to be explored. 

It was argued for a long time that the enhanced value of the critical exponent 
u extracted from the narrowing of the transition region with temperature [34, 122] 
has its origin in the long-ranged disorder present in Ga As-based samples. To our 
knowledge, the present RG approach was the first attempt to quantify this argu- 
ment. As a result the random potential with a power-law correlator leads to values 
of v exceeding ly « 2.35. 

The result indicates that macroscopic inhomogeneities must lead to smaller val- 
ues of K oc l/v. Experimentally, the value of k smaller than 0.42 was reported in a 
number of early (see, e.g., [123] and references therein) as well as recent [124] works. 
This fact was accounted for by different reasons (such as temperature dependence of 
the phase breaking time, incomplete spin resolution, valley degeneracy in Si-based 
MOSFETs, and inhomogeneity of the carrier concentration in GaAs-based struc- 
tures with a wide spacer). Briefly, the spread of the k values was attributed to the 
fact that the temperatures were not low enough to assess the truly critical regime. 
The possibility of having k < 0.4 due to the correlation-induced dependence of the 
effective ly on the phase-breaking length or, ultimately, on the sample size, as in 
Fig. 14, was not considered. 

It should be pointed out that the limited number of RG steps permitted by the 
numerics nevertheless allows one to trace the evolution of the wave functions from 
microscopic scales (of the order of the magnetic length) to macroscopic scales (of 
the order of several /um) which are comparable to the sizes of the samples used in 
the experimental studies of scaling (see e.g., [123, 125]) and much larger than the 
samples [126,127] used for the studies of nicsoscopic fluctuations. Another qualitative 
conclusion of this study is that the spatial scale at which the exponent v assumes 
its "infinite-sample" value is much larger in the presence of the quenched disorder 
than in the uncorrelated case. In fact this scale can be of the order of microns. 
This conclusion can also have serious experimental implications. That is, even if the 
sample size is much larger than this characteristic scale, this scale can still exceed 
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the phase-breaking length, which would mask the true critical behavior at the QH 
transition. 

Let us briefly comment on the possible future application of the method. An ob- 
vious possibility is the application to extended variants of the standard CC model, 
e.g. a network with at least two channels per link in order to describe the mixing of 
Landau levels [53] . There remain several interesting questions such as the behavior 
of the critical properties at the QH transition when changing from strong to weak 
magnetic field. This case could be modelled by bi-directional links in the network, 
which would allow one to trace the transition from the universality class GUE to 
GOE. And last, but certainly not least, we have described the integer QH effect 
within a non-interacting electron picture, but experimental results clearly indicate 
the influence of interactions [31, 128]. Because a full treatment of many-body in- 
teractions is rather difficult one might consider in an approximate view only a few 
interacting particles. In this approach the two-interacting particle problem is re- 
duced to a single-particle problem by increasing the effective spatial dimension and 
including long-range correlations in the disorder potential [129]. Concerning the CC 
model one has to construct an effective four-dimensional network. From this new 
CC network a suitable RG unit should be extracted. This tasks are by no means 
trivial but first steps have already been undertaken [130]. 
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